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ABSTRACT 

Using iV-body simulations of fiat, dark energy dominated cosmologies, we show 
that galaxies around simulated binary systems resembling the Local Group (LG) have 
low peculiar velocities, in good agreement with observational data. We have compared 
results for LG-like systems selected from large, high resolution simulations of three cos- 
mologies: a ACDM model, a AWDM model with a 2 keV warm dark matter candidate 
and a quintessence model (QCDM) with an equation of state parameter w = —0.6. 

The Hubble flow is significant colder around LGs selected in a flat, A dominated 
cosmology than around LGs in open or critical models, showing that a dark energy 
component manifests itself on the scales of nearby galaxies, cooling galaxy peculiar 
motions. Flows in the AWDM and QCDM models are marginally colder than in the 
ACDM one. 

The results of our simulations have been compared to existing data and a new 
data set of 28 nearby galaxies with robust distance measures (Cepheids and Surface 
Brightness Fluctuations). The measured line-of-sight velocity dispersion is oh = 88± 
20 km/sec x (R/7 Mpc). The best agreement with observations is found for LGs 
selected in the ACDM cosmology in environments with — O.i < Sp/p < 0.6 on scales of 
7 Mpc, in agreement with existing observational estimates on the local matter density. 

These results provide new, independent evidence for the presence of dark energy 
on scales of few Mpc, corroborating the evidence gathered from observations of distant 
objects and the early Universe. 

Key words: cosmology: theory - dark matter - large-scale structure of the Universe 
- galaxies: clusters - galaxies: haloes - methods: numerical 



1 INTRODUCTION 

The Hubble flow just outside the Local Group is dynami- 
cally cold, as indicated by the small scatter in the velocity- 
distance diagram of nearby galaxies. This was first pointed 
out by Sandage, Tammann & Hardy (1972). The peculiar 
line-of-sight velocity dispersion <th (see eq. 6 for a defini- 
tion) of galaxies within 5 — 7h~ x Mpc lies in the range 
80 - 40 km s" 1 (Giraud 1986, Schlegel et al. 1994, Sandage 
1999, Ekholm et al. 2001, Karachentsev et al. 2003b). This 
low value is likely linked to the properties of the energy- 
matter density field around the Local Group and it can be 
used to constrain cosmological models. 

Interestingly, cosmological cold dark matter (CDM) 
simulations of hierarchical structure formation have so far 
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not been able to produce values of the velocity dispersion 
around simulated galaxy groups as low as the observed one. 
The three models tested in the works of Schlegel et al. (1994) 
and Governato et al. (1997, hereafter G97) all showed sub- 
stantially larger velocity dispersions and failed to recover 
a significant number of binary systems with the proper- 
ties of our Local Group. Typical values of cth ranged from 
200 km/sec (for the open model) to 500 km/sec (for the old 
"standard" CDM model with Q = 1), several times higher 
than the observational estimates. 

Since these early works, substantial progress has been 
made in the determination of the main cosmological pa- 
rameters. Observations indicate the presence of a significant 
component of smooth energy with large negative pressure 
(observations of high-redshift supernovae by Perlmutter et 
al. 1999, Riess et al. 1999; analysis of fluctuations of the 
cosmic microwave background combined with data on the 
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large-scale structure of galactic distribution by Balbi et al. 
2000, Tegmark, Zaldarriaga & Hamilton 2001, Netterfield et 
al. 2002, Spergel et al. 2003, Verde et al. 2003, Dai, Liang 
& Xu 2004) . This component has been dubbed dark energy, 
and the most appealing candidate is Einstein's cosmological 
constant A (see Peebles & Ratra 2003 for a recent review). 

On theoretical grounds, Baryshev, Chernin & Teeriko- 
rpi (2001) (but see also Chernin 2001) suggested that a 
dark energy (DE) component would cause peculiar veloc- 
ities to cool adiabatically in regions of the universe where 
dark energy overcomes the gravity of the local dark matter 
concentration. In such regions, the usually growing mode 
of the density perturbation is decaying (S g oc a' 1+3lu ' and 
w < —1/3). When dark energy dominates, new structures 
do not condensate and linear perturbations of density and 
peculiar velocities field decay (this effect was also consid- 
ered for vacuum-dominated regions by Chernin et al. 2003) . 
They have shown that for a Milky- Way like halo the DE 
starts to be dominant on a scale of few Mpc and that the 
strength of this effect is related to the equation of state of 
the DE. Their idea was developed further by Karachentsev, 
Chernin & Teerikorpi (2003) and Chernin et al. (2004), us- 
ing a numerical approach to trace the trajectories of the 
radially expanding galaxies back in time in different cosmo- 
logical models to see the effect of a cosmological constant 
on the evolution of the flow. In a more general frame the 
importance of DE on Cluster and Galaxy scales was also 
pointed out by Mota & van de Bruck (2004) and Nunes & 
Mota (2004) using linear theory and spherical collapse model 
and extended to model beyond the standard FRW model by 
Teerikorpi, Gromov, Baryshev (2003). 

Recently, Klypin et al. (2003a) simulated the evolu- 
tion of a cosmic region resembling the real Universe around 
100 Mpc from the Local Group. The density field was 
mapped using the MARKni survey (Willick et al. 1997) and 
a ACDM power spectrum was assumed. The authors report 
a peculiar line-of-sight velocity dispersion of ~ 60 km/sec 
around the Local Group in their simulation, which is a sig- 
nificant improvement compared to previous studies. We have 
extended their work on a bigger sample of simulated LGs to 
test if the low value of an reflects the specific geometry of 
the mass distribution, or/and is a general characteristic of 
the cosmology used. 

The possibility that the observational estimate of an 
adopted in previous works be incorrect or biased cannot be 
dismissed entirely, as the <jh based on Tully-Fisher (TF) dis- 
tances suffers likely from large distance errors (which should 
however increase an) while the lower value estimated by 
Ekholm et al. (2001) was based on a sample of only 9 galax- 
ies. Moreover, the galaxies in the TF sample used by Schlegel 
et al. (1994) had been "regrouped". This was necessary in 
order to avoid inflating an with the internal velocity disper- 
sion of small virialised galaxy associations. No detail, how- 
ever, was given on how the procedure was carried out. If, for 
example, the velocities of loose, non virialised associations 
of galaxies are averaged, the an could be effectively "cooled" 
by removing the peculiar motion of infalling galaxies. 

All considered, the existing literature raises interesting 
questions but leaves them partially unanswered: is the low 
value of an a generic problem of all CDM models or does 
it depend on the cosmology adopted? How important is the 
influence of the local environment and/or the surrounding 



large scale structure? Is the addition of a cosmological con- 
stant important? How robust is the observational estimate 
of an around our own LG? 

Motivated by the growing success of flat A-dominated 
models and the availability of more accurate distance esti- 
mates, we embarked on a new analysis of the problem. In 
this paper, we present a new sample of distances and ve- 
locities of nearby galaxies combined with results from other 
authors on the very local (< 3 Mpc) Hubble flow. We com- 
pare these observational data with new high resolution sim- 
ulations, considering DE models with a constant equation 
of state parameter w different from —1 in order to test the 
impact of DE on the cooling of the local Hubble flow. 



2 REVISITING THE COLDNESS OF THE 
HUBBLE FLOW 

To obtain a robust estimate of the coldness of the Hubble 
flow around the Local Group, it is essential to use accu- 
rate distances to nearby galaxies. As noted above, previous 
measures of an were based on TF distances, which suffer 
from fairly large relative errors, especially for small, poorly 
resolved galaxies and/or viewed at unfavourable angles. 

To this purpose we have gathered data from different 
sources: Cepheid-based distances taken from the Hubble 
Space Telescope (HST) key project to measure the Hubble 
constant (Freedman et al. 2001), distance estimates based on 
the Surface Brightness Fluctuation method (SBF, Tonry et 
al. 2001) and, mainly as a comparison, TF distances (Tully 
et al. 1992). For galaxies common to several samples we have 
kept the Cepheid-based distances. The final sample contains 
28 galaxies within 10 Mpc from the Sun (11 with Cepheid- 
based distances; 17 early type galaxies with distances mea- 
sured using the SBF method: 37 with TF distances). We 
have converted heliocentric distances into distances relative 
to the barycenter of the Local Group, taken 2/3 of the way 
on the line between M 31 and the Milky Way (hereafter de- 
noted as R.2/3)- As for the heliocentric distance to M 31, we 
used the Cepheid derived value of 870 kpc (Ekholm et al. 
2001). The observed heliocentric velocities, taken from the 
NASA Extragalactic Database, were corrected to the values 
they would have if measured from an observer in our galaxy 
at rest relative to the centroid of the Local Group, according 
to Yahil et al. (1977). 

With this sample of 28 galaxies it is also possible to 
measure how cth changes with radius as larger and larger vol- 
umes (with radius up to 10 Mpc) are included. Fig. 1 shows a 
velocity-distance plot, where the distances are relative to the 
barycenter of the LG and the velocities have been corrected 
for the solar motion toward the centroid of the Local Group. 
The crosses show the TF points, the black dots the Cepheid- 
based values and the filled squares are from the SBF sample. 
The straight line corresponds to Ho = 70 km s _1 Mpc -1 . 

Measures of an including the TF sample proved to 
be systematically higher by 30-40% than those using only 
galaxies with Cepheid- or SBF-based distances independent 
of the particular sky quadrant they were selected, showing 
that TF galaxies are not tracing different (hotter) struc- 
tures, but are rather subject to larger distance errors. We 
then decided to not use the TF sample. 

We regrouped galaxies in regions with a number over- 
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density larger than 64. This overdensity threshold corre- 
sponds roughly to virialised region (Eke et al. 1996). Each 
group was assigned the mean position and velocity of its 
members. Only a smaller number (~ 20%) of galaxies were 
linked together, mostly into binary and triple systems. Re- 
sult obtained using our regrouping scheme, do not differ sub- 
stantially from those using galaxy sample without regroup- 
ing. 

The velocity dispersion increases with the size of the 
sphere within which galaxies are included. It varies from 
52 km s _1 within 3 Mpc to 135 km s _1 within 10 Mpc. This 
is in agreement with previous results of Ekholm et al. (2001) 
based on 7 galaxies within 7 Mpc and those of Schlegel et al. 
(1994), based on a sample of galaxies between 3 and 7 Mpc. 
Our estimate of an is a bit higher than the one obtained 
by Karachentsev et al. (2003b) who found <th ~ 42 km/sec 
inside a sphere of 5 Mpc using distances from luminosity of 
the tip of the red giant branch stars of 16 dwarf galaxies, but 
the results are compatible at the 2a level. For an measured 
using HST and SBF galaxies in spheres centered on the LG 
with growing radius R we obtained the following fit: 

a u = 88 ± 20 km/sec x (i?/7Mpc) . 

The final results of the relation oh vs enclosing radius 
are shown in Fig. 5, where we also include the results of 
Karachentsev et al. (2002) for the very local (D < 2.5 Mpc) 
Hubble flow. 



3 A REVIEW OF THEORY PREDICTIONS 

In this section, we very briefly review the analytical work 
done on the subject of cold flows in the context of the Local 
Group Hubble flow. Specifically Baryshev et al. (2001) and 
Axenides & Perivolaropoulos (2002) had already suggested 
on the basis of simple dynamical arguments that dark en- 
ergy could appear as a cooler of peculiar motions in the Lo- 
cal Group environment. What is the expected trend as the 
constant w in the equation of state of dark energy changes? 

A first fundamental difference between the quintessence 
and the A (w = — 1) models is the epoch at which dark 
energy dominates over the gravitational attraction of the 
matter. 

The expansion of a flat universe with dark energy is 
governed by the Friedmann equation 



H(a) = -= H 0V /Q {1 + zf + n DE ,o(l + z) 3 < 1+ "> 



(1) 



where a = 1/(1 + 2) is the scale factor. The dark energy term 
dominates over the gravitational attraction of the matter at 
redshifts lower than z aa : 



1 + z e , 



_ / !2de,o \ 



-l/3w 



(2) 



For w = —1, this gives z cq = 0.326 and for w = —0.6 z cq = 
0.601. 

It has been argued from linear perturbation theory that 
dark energy has little effect on the local dynamics (Lahav et 
al. 1991, Wang & Steinhardt 1998). In linear perturbation 
theory, the peculiar velocity field associated with a spher- 
ical mass concentration with density contrast S averaged 
over proper radius R is t> pcc = where / = fj^f ; 

/ ~ Q a with a ~ 0.6 (Peebles 1980). Wang & Steinhardt 



(1998) have provided an approximation to a in quintessence 
models, which showed that the difference between the linear 
peculiar velocities in models with different w is negligible, 
as the growth factor depends mainly on fio and dark energy 
dominates only at late times. 

However, it has been later recognized that peculiar ve- 
locities can be cooled locally for galaxies located in a "dark 
energy dominated zone" (Baryshev et al. 2001). Consider- 
ing a shell enclosing a mass M evolving in an expanding 
background, one can calculate a "zero-mass" radius beyond 
which the dark energy dominates over the gravitational at- 
traction of the matter (r = 0) . The equation of motion of 
the shell is: 



GM 



±±^V DE ( z )H 2 (z)r. 



(3) 



Baryshev et al. (2001) have argued that peculiar velocities 
are cooled adiabatically for galaxies situated beyond the 
zero-mass surface, where dark energy dominates the dynam- 
ics. For w = — 1, the zero- mass radius is independent of the 
redshift: 



GM 



DE,0-Ho 



= 1.36 Mpc 



(4) 



using M = 2 x 10 12 Mq for the Local Group (van den Bergh 
1999) and our adopted cosmological parameters. For w / 
— 1, we have 



f— 1 

V l + 3w) (1 + z) 



(l + z)(l+«0' 



(5) 



For w = —1/3 (which is equivalent to an open model 
with the same matter energy content), r zm — ► +oo. For 



DE 
zm 



j^ffiirL. This gives rl 



\z = 0) 



w = —0.6, r. 
1.85 Mpc. 

Note that the relative importance of dark energy on the 
dynamics of galaxies around a mass concentration varies in 
time, depending on the local density contrast. The zero-mass 
radius is larger around high-density regions, and for a given 
region it increases with cosmic time. So a galaxy which is 
now close to the Local Group and within the matter dom- 
inated zone could have been in the dark energy dominated 
zone in the past (even if the dynamics of the universe as 
a whole was not dark energy dominated) and have had its 
peculiar velocity cooled during a period of time. Numerical 
simulations make it possible to see the integrated effect of 
the dynamical evolution of galaxies around concentrations 
of mass, which cannot be calculated easily analytically as 
galaxies could move in and out of dark energy dominated 
zones. 

If a galaxy has been in the dark energy dominated zone 
at r > from redshift z cq until now, then its peculiar 

velocity will have cooled adiabatically: 



3(2 = 0) = 



: (^eq) 



1 + z c 



As an upper limit for the peculiar velocity, one can take the 
Hubble expansion velocity vh = HR: 



z (z = 0)<H Q 



2Q Q 



1 + z cq 



flo/lMpc. 



For w = —1, one obtains v pcc (z = 0) < 47.1i?o/lMpc. For 
w = —0.6 one has v pcz (z — 0) < 42.8i?o/lMpc. This upper 
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limit on the peculiar velocity is lower by about 9% in the 
quintessence model with w = —0.6 than in the cosmological 
constant case. 

As pointed out by Baryshev et al. (2001), when w < 
— 1/3 both linear modes are decaying (the second, usually 
growing mode, decays as a 1+ ' iw for z <JC z eq ). Both density 
perturbations and peculiar velocities decay. The lower the 
value of w the faster is the decay, so we would expect to 
have a colder Hubble flow in the ACDM model than in the 
QCDM one, compensating for the earlier dominance of DE 
that we have in models with w > — 1. These two contrary 
effects make hard to predict which Dark Energy model will 
present a colder hubble flow in a given environment. The 
only way to test this scenario is via numerical simulations, 
that we will present in the next section. 

In cosmologies with no dark energy component like the 
open CDM model (equivalent to w = —1/3), the zero-mass 
surface is pushed to infinity and there are no decaying modes 
so peculiar velocities can't be cooled, which is consistent 
with the high peculiar velocities seen in simulation of the 
OCDM model by Governato et al. (1997). 

Finally, Axenides & Perivolaropoulos (2002) have calcu- 
lated the growth suppression factor of velocity fluctuations 
in a dark energy model compared to the Einstein-de Sitter 
model and concluded that the presence of dark energy is not 
sufficient to explain the coldness of the local flow. However, 
this holds true for cosmological models with different w and 
the same high redshift normalization of the density field per- 
turbations. However, realistic models must be normalized to 
match the observed value of cr 8 at z = 0. 

From the above analysis it seems plausible that cos- 
mologies with a dark energy component would present colder 
flows around cosmic haloes. However the magnitude of this 
effect, the dependence on w and on the local properties of 
the density field remain poorly constrained by simple ana- 
lytical arguments based on linear extrapolations. These con- 
siderations make resorting to numerical iV-Body simulations 
necessary. 



4 SIMULATIONS 

We present results from four simulations, performed using 
two different codes: 

- ACDM: two simulations, both with h = 0.7, Q = 0.3, 
Q,\=Q.7, as — 1.0 but with different resolutions (see details 
below) ; 

- AWDM: a flat model similar to the low resolution 
ACDM but with a sterile neutrino mass of 2 keV as a warm 
dark matter candidate; 

- QCDM: a quintessence model with the same param- 
eter of the previous two simulations AC(W)DM model but 
with a constant equation of state parameter for the dark 
energy w = Pq/pq = -0.6. 

The A simulations were performed using PKDGRAV, a par- 
allel, multistepping tree-code with periodic boundary condi- 
tions (Stadcl 2001). For the quintessence simulation, we used 
the parallel Adaptive Refinement Tree code ART (courtesy of 
A. Klypin, Kravtsov, Klypin & Khokhlov 1997). 

QCDM models with w > —0.7 are now disfavored by 
CMB observation (Kogut et al. 2003, Verde et al. 2003). 
Nevertheless we have used this value for w in order to try 



to maximize the influence of dark energy on an (so it will 
be used as a limiting case). According to theoretical expec- 
tation (Baryshev et al. (2001) and section 3) the strength of 
this effect depends on the epoch when the DE becomes the 
dominant component of the universe: the earlier the DE pe- 
riod begins, the stronger the cooling. In our QCDM model 
the DE period starts earlier than in the ACDM model, and 
so we would expect a larger effect of dark energy (but see dis- 
cussion in earlier paragraph). To calculate the power spec- 
trum P(k) for the AWDM model and the QCDM used the 
CMBFAST code to generate transfer functions (Seljak & Zal- 
darriaga, 1996). Sterile neutrinos are a popular candidate 
for warm dark matter (Dodelson & Widrow 1994, Dolgov 
& Hansen 2002). Though formation mechanisms for sterile 
neutrinos are likely to produce a non-thermal energy spec- 
trum, the filtering effect may still be modeled accurately 
using thermal neutrinos with adjusted masses (Hansen et 
al. 2001). The AWDM model can be characterized in terms 
of the reduction in fluctuations on dwarf galaxy scales. For 
2 keV neutrinos the fluctuations on the scale of 10 11 Mq 
(00.59) are down by 10% over ACDM. The power spectrum 
at 2iv/k — r = 0.59fe _1 Mpc is lower by 2 orders of magni- 
tude and thus contribution of small objects to U0.59 is vastly 
reduced. Fluctuations on the scale of 10 12 Mq haloes are 
only marginally affected (<7i.26i4i down by 3%). 

Each AC(W)DM simulation cube is 100 Mpc on a side 
with both using 144 3 particles and a spline kernel softening 
of 60 kpc (the spline kernel is completely Newtonian at 2 
softening lengths) The QCDM simulation has a box size of 
90 Mpc with 128 3 particles, which gives roughly the same 
mass resolution as the low resolution A run. Since our goal 
is to identify galactic sized dark haloes, not to resolve their 
internal structure, this combination of mass/force resolution 
is adequate. Each simulation used at least 512 steps for ev- 
ery particle, and many more (up to a few thousands) for 
particles with large accelerations as those at the center of 
dense haloes. The mass resolution is such that a well resolved 
halo (30 particles) at redshift zero has a circular velocity of 
95 km/s. We can therefore resolve isolated haloes that are 
associated with galaxies that have luminosities as low as a 
few % of an L* galaxy (Cole et al. 2001). 

We have also analysed a very high resolution simula- 
tion for the ACDM model, with 81 x 10 6 (432 3 ) dark matter 
particles simulated from a starting redshift zo = 69. The 
employed volume is 50 /i _1 Mpc, which gives a particle mass 
of 1.3 x 10 8 /i _1 Mq allowing us to resolve haloes down to less 
than 10 10 /i _1 Mq with 75 particles. The force resolution is 
5 IT 1 kpc and the particles in the highest density region un- 
dergo 20 000 timesteps (this simulation is described in more 
details in Reed et al. 2003a). In this run, we are able to re- 
solve the substructure inside haloes corresponding to galaxy 
groups, allowing us to evaluate the effect of substructure 
within larger haloes in velocity dispersion measurements. 
The results of a comparison with the low resolution run are 
shown in Fig. 3. 

Our ACDM model has parameters very close to those 
currently suggested by observational results. This model has 
been successful on an extended number of crucial tests on 
scales down to a few Mpc. However, recent results pointed 
out a possible flaw of ACDM at galactic and subgalactic 
scales, like steeper density profiles than those implied by 
galaxy rotation curves (Moore et al. 1998, de Block & Bosnia 
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2002, Swaters et al. 2003, Reed et al. 2003b, Diemand et al. 
2004a) and an excess of substructure (i.e. small satellites) in 
galactic haloes (Moore et al. 1999, Klypin et al. 1999, Ghigna 
et al. 2000, De Lucia et al. 2004, Diemand et al. 2004b; 
Klypin et al. (2003b) for QCDM model). Both problems may 
be linked to an excess of power at small scales. Hence, as a 
cure to the above mentioned problem a WDM rather than 
CDM component reduces the amount of density fluctuation 
and subgalactic scale and may prove useful in cooling the 
Hubble flow around the Local Group. To boost its possible 
effects we used for the mass of the WDM a value of 2 keV, 
this is close to the minimum allowed from constraints based 
on satellites abundance and the abundance of absorption 
systems (Bode, Ostriker & Turok 2000) . A lower WDM mass 
would likely underproduce both of them. Recently Yoshida 
et al. (2003) put strong constraints on WDM candidates 
lighter than lOkeV based on the abundance of haloes at the 
epoch of reionization. 
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5 HALO FINDERS AND SELECTION OF LG 
SYSTEMS 

Our simulations only follow the evolution of the dark matter 
component. It is expected that baryons condense and form 
galaxies at the centers of the dark matter haloes (White & 
Rees 1978; Cen & Ostriker 1992,1993a,b; Katz, Hernquist & 
Weinberg 1992). For our purpose of locating LG halo can- 
didates (not resolving their internal structure), we resorted 
to two algorithms: SO and SKID (SKID is publicly available 
at http://www-hpcc.astro.washington.edu). Our version of 
SO (Lacey & Cole 1995) uses SKID halo centers (based on 
the density maxima) as a start, then a sphere is grown 
around this center until the enclosed density drops below 
some threshold which is chosen to be that associated with 
the average density inside a spherical virialised halo in that 
given cosmology (values were 97.1p cr it and 135.2p cr it for 
LC(W)DM and QCDM respectively (Mainini et al. 2003); 
At difference with SKID, SO imposes a spherical symmetry 
to the haloes and, perhaps more importantly, regroups sub- 
structure within larger virialized structures. We verify that 
this does not affect our results in subsection 6.2. We veri- 
fied that our results were largely independent on which halo 
finder we used and decided to show results only from SO. In 
this paper we assume that a galaxy with a certain circular 
velocity would be found in a dark halo that has a similar 
circular velocity (Cole et al. 2001). 



Figure 1. Expansion velocity Vlg in the LG reference frame as 
a function of radial distance Rlq from the Local Group baricen- 
ter. Black dots: Cepheid based distances; Filled squares: Surface 
Brightness Fluctuations based distances; Crosses: TF distances. 
The straight line is for Hq=70 km/sec/Mpc. 



(ii) A LG sample defined such that s < 1.0 Mpc, circular 
velocities 125 < V c < 270 km/sec, negative radial velocities 
and no neighbors within 3 Mpc with circular velocity larger 
than either of the two group members. 

(iii) A LG sample defined as (2) but with the addi- 
tional requirement that the binary haloes must lie close 
(5-12ft" 1 Mpc) to a Virgo sized cluster with 500 < V c < 
1500 km/sec. 

In the ACDM cosmology, about 30 Local Groups satisfy 
criteria 1) and 2), and about 15 satisfy point 3) as well. 
As we will show later there are no significant differences 
between samples 2 and 3 (and actually with the global field 
halo population, see Fig. 6), and we will refer to our Local 
Group sample as systems identified using criterion 2). 

We define the peculiar line-of-sight velocity dispersion 
from the Hubble flow for a set of haloes within a sphere of 
radius R (in Mpc) centered on the center of mass of each 
simulated LG as: 



6 STATISTICAL PROPERTIES OF LOCAL 
GROUP CANDIDATES: COMPARING 
SIMULATIONS WITH THEORY AND 
OBSERVATIONS 

6.1 Selection of Local Groups 

Following G97, from the halo catalogs obtained we selected 
"Local Groups", i.e. binary haloes with the following crite- 
ria: 

(i) A generic sample of binary haloes with separations 
s < 1.5 Mpc and circular velocities 125 < V c < 270 km/sec. 



oh = 



n 

\ n i=i 



where 



71 

(v - H loc D) = i - HiocA) 



(6) 



(T) 



where a local Hubble constant Hi oc is fit to each dataset of 
haloes. We have used the same procedure also to analyse the 
data points shown in Fig. 5. 
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Figure 2. Hubble flow velocity dispersion in a sphere of 10 Mpc: 
ACDM (circles), A WDM (triangles) and QCDM (squares). The 
lines marks the results from models described in G97, the large 
starred dot marks the real Local Group. 



6.2 Comparison with observations 

Comparing the results from this new set of simulations to 
the ones obtained by G97 (see Fig. 2), it is evident that the 
new, dark energy dominated models fare better, expecially 
when compared to the new observational estimate of an- 
The velocity dispersion is down by at least a factor of two 
compared to the classical, but outdated SCDM model for a 
large range of overdensities. 

ACDM and AWDM give almost the same oh for the 
whole range of overdensities. The QCDM model produces a 
colder Hubble flow, as better shown in Fig. 4. 

In Fig. 3 we show the comparison between low and high 
resolution simulations. As expected, the high resolution run 
has a higher value of an at small scales, due to the fact 
that we are able to resolve the inner structure of the more 
massive haloes (corresponding to galaxy groups) around the 
LGs. The velocity dispersion is lower on scales > 7 Mpc. 
This could be related to the smaller box size that reduces the 
effect of the large-scale structure. Nevertheless, both simula- 
tions agree within the la error bars showing that the effect 
of substructure is small. 

As shown by G97, the velocity dispersion around the 
Hubble flow depends strongly on the local overdensity. This 
can be determined for the Local Group using the complete 
catalogue of IRAS galaxies. The current best estimate based 
on IRAS galaxies is Sp/p — 0.60 ± 0.15 for a top-hat sphere 
of radius 500 km/s (Strauss 1996, private communication). 
Hudson (1993) uses a compilation of optical galaxy surveys 
to study the local density field and within the same volume 
he finds an overdensity of ~ 0.2. The agreement between 
the optical and IRAS surveys is encouraging and indicates 
that the Local Group resides in a region of the Universe of 
almost average density. 



Figure 3. Hubble flow velocity dispersion computed in spheres 
of growing radius for the high-resolution ACDM simulation. 




4 6 8 10 




Sphere Radius (Mpc) 

Figure 4. Hubble flow velocity dispersion computed in spheres of 
growing radius for different cosmological models. The lower panel 
shows results for LGs with —0.1 < Sp/p < 0.6, the upper panel 
for all the field galaxy-sized haloes. 

If simulated LGs are selected in regions with a local 
overdensity between —0.1 < Sp/p < 0.6, both the normal- 
ization and the shape of the an vs R relation are well re- 
produced (see Fig. 5) - a remarkable success for flat, A- 
dominated models. 

The AWDM performs equally well (cf Fig. 4, the higher 
value of an at very small scales is of small statistical signif- 
icance), indicating that it is the presence of a cosmologi- 
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Figure 5. Hubble flow velocity dispersion computed in spheres of 
growing radius. The large filled dots mark results from our new 
estimate of oh f° r the Local Group, the last starred point was 
taken from the estimate by Karachentsev et al. (2002). Simulation 
results are shown for the ACDM model for different choices of the 
density environment. Note that the flattening of the dotted curve 
at small Rs is due to the low number of LGs in that overdensity 
range. The first two points (R=2 and R=3) are made both with 
just 2 LGs, so the statistical variance is quite high. 



cal constant, which is the most important ingredient, rather 
than the details of the power spectrum at subgalactic scales. 
This is further supported by the QCDM simulation, which 
produces a Hubble flow colder than in both A models at all 
scales. That seems to indicate that a longer period of dark 
energy dominance has the capability to produce a colder 
Hubble flow, in agreement with the theoretical expectations 
by Baryshev et al. (2001) and section 3). Comparing our 
analysis to the results of G97, it appears that only DE- 
dominated models are able to satisfy at the same time small 
and large scale structure constraints and generate a sub- 
stantial number of simulated LGs with <th comparable to 
the observed values. Also, LGs in LC(W)DM cosmologies 
fail to have <th in the correct range if their local overdensity 
is too high or too low, as shown in Fig. 5. 

How are important the details of the mass distribution 
around the LG other than the average density on a few Mpc 
scale? Klypin et al. (2003a) obtained a cold Hubble flow 
around a simulated region that replicated the mass distri- 
bution inferred from a redshift survey of the local Universe. 
Their simulation assumed a ACDM cosmology. Their work 
was based on one constrained simulation of our local en- 
vironment, our LG samples were selected using some sim- 
ple selection environment criteria like the absence of nearby 
massive groups or the presence of a rich cluster within 5- 
12 Mpc. However, no substantial difference was found be- 
tween the typical <th in the different subsamples depending 
on the presence of nearby clusters. 

To further test the hypothesis that a special local mass 



Figure 6. Peculiar velocity around LG candidates and field 
galaxy-sized haloes (small dots) in the ACDM model 



distribution is necessary for a cold Hubble flow, we repeated 
our analysis for a much larger sample of Milky Way (MW) 
sized field haloes, with no special selection criteria. Results 
are shown in Fig. 6. A generic sample of MW-sized field 
haloes shows values of <th in the correct range if the local 
environment density is around average. This shows that the 
Local Group environment and details of its mass distribu- 
tion are not special in any significant way and that in DE- 
dominated cosmologies a small oh is a general characteristic 
of the environment around field galaxies. 



7 CONCLUSIONS 

Observational evidence at large and intermediate scale 
points to a flat, A-dominated CDM Universe. The coldness 
of the Hubble flow provides another, independent test at in- 
termediate scales, where different kinds of dark matter and 
energy could play a significant role in shaping the local dy- 
namics. 

In this paper we have tested the predictions from cos- 
mological models now strongly supported by observational 
evidence against a new, robust estimate of the Local Group 
<th in spheres of increasing radius up to 10 Mpc based on a 
sample of nearby galaxies. 

We have shown that simulations are able to successfully 
reproduce the observed cold Hubble flow around field galax- 
ies only in a flat, DE (w < —1) dominated cosmology and 
only in environments where the local density measured on a 
scale of a few Mpc is close to the average one. The inclusion 
of a warm dark matter component does not change this re- 
sult significantly, as removal of small-scale power does not 
play a major role. 

We have also shown that w > — 1 (quintessence mod- 
els) where the era dominated by dark energy begins earlier, 
produce an Hubble flow marginally colder than ACDM. 
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This is in agreement with the theoretical expectation 
(see section 3 and also Baryshev et al. 2001, Chernin et 
al. 2004) that the coldness of the local Hubble flow follows 
from the fact that a galaxy's peculiar velocity decreases via 
adiabatic cooling if it spends a sufficient time in the dark 
energy dominated region where linear perturbations decay. 
The presence of a dark energy field is then the key ingredient 
to explain the observed coldness of the local Hubble flow. 

While the right cosmology and environmental density 
are crucial, the local dynamical criteria do not single out the 
LG from the generic population of field galaxies. A cold Hub- 
ble flow around Milky Way sized galaxies is equally common 
as in our LG subsample. This implies that the details of the 
mass distribution in our local environment, like the vicinity 
to larger clusters and/or voids and filamentary structures, 
do not play a major role in shaping its velocity field. 

The environment of our Local Group provides a new, 
independent evidence for the existence of dark energy on 
scales of just few Mpcs, corroborating the evidence gathered 
from observations of distant objects and the early Universe. 
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